SU(2) instantons with boundary jumps and spin tunneling in magnetic molecules 
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Coherent state path integrals are shown in general to contain instantons with jumps at the bound- 
aries, i.e., with boundary points lying outside classical parameter or phase space. As an example, 
the magnetic molecule Fes is studied using a realistic Hamiltonian, and instanons with jumps are 
shown to dominate beyond a certain external magnetic field. An approximate formula is found for 
the fields where ground state tunneling is quenched in this molecule. 
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The purpose of this paper is to discuss two problems, 
a specific one and a general one. The specific problem 
concerns the tunneling between ground Zeeman levels of 
certain magnetic molecules Q . The general problem con- 
cerns the nature of tunneling paths, or instantons, in co- 
herent state or phase space path integrals [0 . As a rule, 
such paths can be found only if one complexifies phase 
space, i.e., allows either the momenta, or the coordinates, 
or both, to become complex. However, in all previous 
studies of which we are aware, the instantons always start 
and end in real phase space, at the points corresponding 
to the classical energy minima. The new instantons re- 
ported here, by contrast, do not even have end points in 
real phase space. We refer to these as boundary jump 
instantons. Such paths are a basic part of the formal 
structure of coherent state and phase space path inte- 
grals H U, but there has never been a need to include 
them in tunneling problems, as paths without jumps have 
always been available ||. This, as we shall show, is not 
an accident. Boundary jump instantons are analogous to 
extra saddle points in the method of steepest descents 
for one-dimensional integrals, and like them, may or may 
not be relevant in any given situation. But, one can not 
ignore them a priori. In this paper we give the rules for 
finding these extra paths and their contribution to tun- 
neling, and apply them to our illustrative example. 

Our work is part of the broader program of studying 
the semiclassical limit of spin systems. Tunneling is just 
one of the phenomena amenable to semiclassical meth- 
ods. Many problems, from rotating molecules to many- 
body aspects of nuclear structure, can be modelled in 
terms of a large spin, for which the semiclassical approx- 
imation is the natural one 01. It is perhaps not widely 
realized that this limit is not as well understood for spin 
as it is for massive particles. The correct semiclassical 
spin propagator has been recognized only gradually since 
the late 1980's [0, and its consistency under composition 
of successive propagators was shown only recently || . 

We begin by describing our specific problem. The 
molecular ion [(tacn) 6 Fe 8 02(OH)i 2 ] 8+ (or just Fes for 
short) forms a solid in which the Fes groups are essen- 
tially noninteracting, and each behaves like a single spin 
of magnitude J = 10 in its ground manifold. The degen- 



FIG. 1: Re S for the Fes instantons, marked 1-4 as shown. 
S3 and S4, are purely real. Also shown is the ground pair 
splitting A(h). 



eracy of the 21 Zeeman levels is partly lifted by spin-orbit 
effects. Including an external magnetic field, the system 
is well described by the anisotropy Hamiltonian |9j 



n = kxJt 



C[Jl + J^]-gii B J z H, (1) 



where J = (J x , J y , J z ) is a spin operator, and J± = 
J z ±iJ y (not J x ±iJ y ). If H = 0, the spin has degenerate 
classical minima along the ±x axes, which cant symmet- 
rically toward z as H is turned on. We wish to study 
the splitting A due to quantum tunneling between these 
minima, particularly as a function of H . It is convenient 
to work in reduced variables A = k 2 /ki, A 2 — CJ 2 /ki, 
and h = H/H Cl with H c = 2k\J / g[iB- Measurements 
yield g ~ 2, k x ~ 0.338 K, k 2 ~ 0.246 K, and C ~ 29 

The splitting A (H) oscillates with H, being quenched 
at certain H p], II . The effect can be viewed in terms 
of two instantons that wind about the hard axis, z, in 
opposite directions, and interfere with one another [ p^[ . 
Their actions differ by i J A, where A is a real Berry phase 
equal to the area of the closed loop formed by the two 
paths on the complexified unit sphere. A(H) = when- 
ever JA(H) is an odd multiple of ir. If C = 0, there 
are 2 J quenching fields, at H = ( J — n — 1/2) AH, with 
n = 0, 1, . . . , 2J-1, and AH = (1-X)^ 2 H C /J = 0.263 T. 
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Turning on even a small C has a big quantitative effect. 
With the value for Fes, e.g., direct numerical diagonaliza- 
tion of TL (see Fig. [l]) reveals only 4 quenching points for 
H > instead of 10. The spacing remains quite regular, 
but is increased by ~ 50%. Given the visual appeal of 
the interfering instanton picture, it is interesting to ask 
how this picture is modified when C =/= 0. To this end, 
we review some aspects of SU(2) instantons. 

In the instanton method, one seeks a propagator such 
as Kfi — (zf\ exp[— TIT] \zi) via the path integral 



K 



fi 



d[z]d[z]e- s W>>*V» 



(2) 



in the limit T — > oo. Here, \z% f) are (unnormalized) 
spin coherent states defined for any complex number z by 
\z) = e zJ - | J, J), where | J, J) is the eigenstate of J 2 and 
J z with eigenvalues J(J + 1) and J. Further, J • n\z) = 
J\z), if z = tan(6>/2)e l( ^, (8, </>) being the spherical polar 
coordinates of n. We will take the points z, and z/ to 
be degenerate minima of the classical energy. S is the 
action for a path specified by z(t) and z(t), and is given 
by 



S=- 



T/2 r 



T/2 



J- 



1 



H(z,z) 



dt. 



(3) 



Here, z and z are formal complex conjugates, but should 
be regarded as independent variables, since both are 
needed to determine 8 and <fr, e.g. Further, 



H(z',z) = {z'\H\z)/{z']z). 



(4) 



For TVs such as (|T|) that are polynomial in Jj, H(z',z) 
is holomorphic in z and antiholomorphic in z' ]l4| ]. In 
Eq. (||), the first term is the Wess-Zumino or Berry phase 
term. We will refer to the two terms in S as the kinetic 
and dynamical terms, Sk and Sn- 

Instantons are paths that start at Zj and end at Zf, 
and obey the Euler-Lagrange (EL) equations, 



(1 



dH 



(1 



dH 



2 J 



dz 



2 J 



dz 



(5) 



It is easily verified that along these paths energy is con- 
served, i.e., dH(z,z)/dt = 0. Because of this, and be- 
cause Zi and Zf are energy minima, one cannot find a 
solution lying on the real unit sphere. We must allow z 
and z to be completely independent complex variables. 
In other words, z(t) need not equal z*(t), where the 
star denotes the true complex conjugate. However, one 
can always find an instanton with endpoints on the real 
sphere. This is because, if we denote the minimum en- 
ergy by E mm , the instanton obeys H(z(t), z{t)) = E m i n , 
and one clearly has H(z*, Zj) = H(zp Zf) — E m [ n . 
Let us illustrate this using Eq. p) with C = 0. Then, 



H(z, z) = fciJ 2 



(1 - zzf - A(z - zf - 2h(l - z 2 z 2 ) 



(1 



(6) 



The minima are at z = z = ±zq where zo = [(1 — + 
h)] 1 / 2 , and E min — —k\J 2 h 2 . From H(z,z) = E min we 
obtain z(z): 



VAz ± (1 - h) 
VA± {l + h)z 



(7) 



These equations give the instanton trajectories in z-z 
space, without giving the t dependence. For general z, 
z z* , but if z = ±zo, z = z* . Thus the instanton 
endpoints are on the real sphere. From Eq. (^), we can 
now evaluate S, and recover previous results, along with 
the interference effect |TT| . 

If we now turn on C, the solutions (^) will evolve 
smoothly, and continue to have classical end points. They 
will continue to interfere, and one can find the fields 
where A vanishes by calculating J A numerically. When 
we do this, wc find that the spacing agrees closely with 
the answer from direct diagonalization oiTL, but we also 
find, incorrectly, quenching points at h > ho ~ 0.25. 

The problem is that we have not formulated the prin- 
cipal of least action (or Hamilton principal function, to 
be precise) sufficiently carefully H |^, |). One must in 
fact include an explicit boundary term Sb in S: 



Sb = Jin 



(l + z(-T/2) Zi )(l + z f z(T/2)) 

(l + Z*Zi)(l + Z f Z* f ) 



(8) 



If we now vary S = Sk + So + Sb including the end- 
points, and set SS to 0, we discover of course the EL 
equations ([5]), but also that SS has no terms in 5z{— T/2) 
and 5z(Tj2). This means that the boundary conditions 
on Eq. (F 



z{-T/2) = z h z{T/2) = z f , 



(9) 



and that z, = z(—T/2) and Zf = z(T/2) must be left free. 
Otherwise, we would have four boundary conditions on 
a second order system of differential equations, and the 
problem would be overdetermined. The term Sb can 
also be found by careful time-slicing of the propagator. 
Its inclusion in S has many other nice consequences: e.g., 
the Hamilton- Jacobi equations 



dS cl 



= 2.J- 



z f 



1 



dS c 

dz. 



= 2J- 



1 



(10) 



Since Zi and zt are not determined, one may have so- 
lutions to Eqs. (^) and (||) with z\ ^ z*, Zf ^ Zj. These 
are the boundary jump instantons. Their velocities z and 
z do not vanish at the end points because the derivatives 
dH/dz and dH/dz are not zero. Thus the instanton du- 
ration is finite, and although the energy E — H{z, z) is 
still a constant of motion, the value of E is not imme- 
diately obvious. In fact, since Sd = J Hdt, wc must 
choose E = -Emin- Otherwise, when we sum multiin- 
stanton terms, the instantons with jumps will trivially 
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dominate or be dominated by those without jumps. This 
point comes out more easily in Klauder's formulation. 
He argues that since the continuum path integral is a 
formal construct with meaning only as a limit of its dis- 
crete version, one may add a term to the integrand for 
S that is quadratic in the velocities z and z, with an in- 
finitesimal coefficient e that is sent to at the end. The 
EL equations are then a fourth order system, and one 
may specify all four zi, z~i, Zf and zj. The instantons 
with jumps then appear as solutions to the EL equations 
with internal boundary layers of thickness 0(e) since the 
terms in z and z have coefficients e. Energy is conserved 
in these boundary layers too, and when one takes the 
e — > limit, they yield a contribution that is explicitly 
independent of e (which makes the procedure legitimate), 
and is precisely equal to Sb above. Note that Sb = for 
an instanton without jumps. 

Hence the general procedure for finding all instantons 
is as follows. For any Ti. with degenerate minima, we first 
find i? m in, and the classical minima (z*, 2,), {z**, zf). We 
then find the allowed values of Zi by solving 



H(z, z^ = E min . 
This equation has a double root at z = z*, since 



d_ 

dz 



H(z, z) 



0_ 

dz 



H(z,z) 



= 0. 



(11) 



(12) 



However, it may also have additional roots at z ^ z* , 
which will then be the end points of instantons with 
boundary jumps. (A completely analogous procedure ap- 
plies to Zf.) We then obtain z(z) for all instantons from 
energy conservation, making sure that they connect on to 
the appropriate end points. This is enough to compute 
Sk and Sb for each instanton (the time dependence is 
not needed) , and Sd — E m - m T for all of them. If we label 
the various instantons by a, we can write 
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FIG. 2: A jump instanton (number 3), for H = 0.1H C . Suf- 
fixes e, m, and h denote easy, medium, and hard axes. 



z(z) has four branches, correpsonding to the different in- 
stantons. We label the first two instantons, which have 
Im S ^ 0, and interfere with each other, 1 and 2, and 
the last two, which have jumps and Im S — 0, 3 and 4. 
An instanton with a jump at one end also has a jump 
at the other end. A 180° rotation about z sends instan- 
ton 1 into 2, which guarantees Re Si = ReS 2 , 71 = 72- 
We show instanton 3 in Fig. |2|, 1 and 2 in Fig. || (all for 
h = 0.1), and ReS a (h) in Fig. |l|. For any h, the domi- 
nant instanton is that with the least Re S. Hence, except 
in the immediate vicinity of ho, only instantons 1 and 2 
are relevant for h < ho, and only 3 is relevant for h > ho . 
This explains why A does not oscillate for h > ho- 

We can find the quenching fields numerically, but we 
have also found an analytic approximation, based on the 
small parameter £ = 4A2/1 2 , which explains why they are 
so regularly spaced. This result may also be of wider 
interest, since A oscillations have now been seen in an- 
other system [M. To derive this, it is better to use polar 
coordinates. With u = cos 9 and s = sin0, the energy 



A = J2 



7c e 



(13) 



where 7 a is the prefactor arising from integrating over 
Gaussian fluctations about each instanton. On physical 
grounds we expect 7 Q to be of the same order for all a 
for smooth Hamiltonians, and it may be estimated as the 
small oscillation frequency about the minimum. (For in- 
stantons without boundary jumps, Ref. pl| formulates 
how to find j a .) Hence, the relative importance of vari- 
ous instantons is determined largely by the actions S a . 

Let us now return to our model (0) . It may be verified 
that H(z,z) = P(z,z)/(\ + zz) A , where P is a polyno- 
mial of degree 4 in z and also in z. Thus Eq. ( pi] ) is a 
quartic in z. Two of its roots are indeed z* , and connect 
on to instantons without jumps, but two are different 
and distinct and connect to instantons with jumps. The 
equation H(z, z) = E m { n is also a quartic and the solution 




FIG. 3: Instantons without jumps for H = Q.1H C . We show 
\Ji\, since now all Ji are complex. All are even about 
t = 0, and Jm slowly asymptotes to as t — » ±00. 
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conservation condition takes the form 

g(u, s) = -\Z{s)u A + R(s)u 2 - 2hu + W(s) = 0, (14) 
where 



Z(s) 
R(s) 
Wis) 



4A 2 (l + 6s 2 + ,s 4 ) 



1 - As 2 



12A 2 s 2 +4A 2 s 4 , 



go + h 2 + Xs 2 - 2A 2 s 4 , 



(15) 
(16) 
(17) 



with <7o = — (A + h 2 ) — E mul sa 2A 2 /i 4 . g(u,s) has four 
roots u(s). Regarding s as real, we are interested in the 
complex conjugate pair of roots which tend to the energy 
minima 6 = 9 , as s — > 0. Let the real and imaginary 
parts of these roots be A(s) and B(s), i.e., let 

u(s)=A(s)±iB(s). (18) 

From the imaginary part of Eq. ( |l4| ) we get 

B 2 = A 2 - (R/Z) + (h/AZ), (19) 

and if we substitute this result for B into the real part, 
we get an equation for A alone: 



4Z 2 A 6 - ARZA 4 + (R 2 + 2WZ)A 2 - h 2 



0. 



(20) 



We now make the self-consistently verifiable assumption 
that A = 0(h). Then the terms ZA 4 and Z 2 A 6 are 
O(C) and 0(( 2 ) relative to the the remaining terms, and 
may be dropped. This yields A = h(R 2 + 2WZ)~ 1 ^ 2 . 
The quantity R 2 + 2WZ can be seen to be a fourth or- 
der polynomial in s, and depends on h only through the 
combination A 2 (/i 2 + go), which is 0((). If we neglect 
this weak h dependence, we get 



A(cf>) 
Pi 
Pi 



h(l + P 2 sin 2 </> + P 4 sin 4 </> 

-2A + 24A 2 + 8A 2 A, 

A 2 + 8A 2 + 24AA 2 + 128A 2 



-1/2 



(21) 
(22) 
(23) 



Since Sk — ^/(l 

r 2TT 



cost 



dt in 9, 4> variables, A = 
2tt — Jq V A(4>)d<p. If we keep only instantons 1 and 2, 
A = 271 exp(— Re5*i) cos( J A/2) up to a phase factor. 
Using Eq. (pl|), we see that the zeros of A are equally 
spaced with spacing AH = ttH c / JI(X, A 2 ), where 



/(A,A 2 



(1 + P 2 sin 2 <j) + Pi sin 4 



a/2 



(24) 



For the Fes parameters, / = 3.88, implying AH = 
0.409 T. The experimental value is 0.41 T. 

We conclude with some general remarks about instan- 
tons with boundary jumps (or internal boundary layers, 
in the Klauderian view). It is clear that they must be 
present in all coherent state path integrals, not just for 
spin, and our discussion is easily extended to these cases. 
It would be interesting to find other concrete instances 



where they occur, both in quantum mechanics, and in 
field theories. It would also be interesting to reaxamine 
problems such as a particle in a one dimensional potential 
well in a coherent state formulation. 
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